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I. INTRODUCTION 



The search for one or more Higgs bosons is the central task of the CERN Large Hadron 
Collider (LHC), with y/S = 14 TeV and a luminosity of 100 fb" 1 per year. In the Standard 
Model (SM), the Higgs boson mass is a free parameter with an upper bound of mj? < 600 
- 800 GeV Beyond the SM, the Minimal Supersymmetric Standard Model (MSSM), 
whose Higgs sector is a special case of the Two Higgs Doublet Model (2HDM) 3], is of 
particular theoretical interest, and contains five physical Higgs bosons: two neutral CP-even 
bosons h° and H°, one neutral CP-odd boson A , and two charged bosons H ± . The h° is 
the lightest, with a mass m h o < 140 GeV when including the radiative corrections Q], and is 
a SM-like Higgs boson especially in the decoupling region (m^o ^> m^o). The other four are 
non-SM-like ones, and the discovery of them may give the direct evidence of the MSSM. It 
has been shown in , |5| that the h° boson of MSSM cannot escape detection at the CERN 
Large Hadron Collider (LHC) and that more than one neutral Higgs particle can be found 
in large area of the supersymmetry (SUSY) parameter space 

At the LHC, the neutral Higgs bosons can be produced through following mechanisms: 
gluon fusion gg — > B fl 33, weak boson fusion qq - qqV*V* -> qqW/qqlP Q, 
associated production with weak bosons associated production with a heavy 

quark- ant iquark pair gg,qq — > ti<fi/bb(f) jlJ] and pairs production Studying the asso- 

ciated production process of a neutral Higgs boson and a vector boson at future hadron 
colliders may be an interesting way in searching for neutral Higgs bosons, since the total 
cross section may be large and also the leptonic decay of the vector boson can be used as a 
spectacular event trigger. In the SM, the process qq^ — > W/Z°h° SM has been studied both 



at the leading order (LO) |ll( and the next-to-leading order (NLO) 13|, |l6j in QCD. In 



the 2DHM and MSSM, the associated production of h°(H°)Z° and A°Z° has been studied 



only at tree level for Drell-Yan process and at one-loop level for gluon fusion in [12 1 and 
H H Q , respectively. 

It was shown in Ref. that the A°Z° associated production rate at the LHC strongly 
depends on the SUSY parameters tan (3 (the ratio of two vacuum expectation values) and 
rriA (the mass of A ). The total cross section increases with increment of tan/3, and de- 
creases with increment of tua- In this paper, we present the complete NLO QCD, including 
supersymmetric QCD, calculation for the cross section of the associated production of A°Z° 
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through bb annihilation process at the LHC For simplicity, in our calculation, we neglect 
the bottom quark mass except in the Yukawa couplings. Such approximations are valid in 
all diagrams, in which the bottom quark appears as an initial state parton, according to the 
simplified Aivazis-Collins-Olness-Tung (ACOT) scheme [21]. To regularize the ultraviolet 
(UV), soft and collinear divergences, two regularization schemes are used in our calcula- 
tions for cross check, i.e. the dimensional regularization (DREG) scheme 221 ] (with naive 
75 scheme (23 1) and the dimensional reduction (DRED) scheme 24], and their results are 
compared. 

The arrangement of this paper is as follows. In Sect. II, we show the LO results and 
define the notations. In Sect. Ill, we present the details of the calculations of both the 
virtual and real parts of the NLO QCD corrections, and compare the results in DREG with 
those in DRED. In Sect. IV, by a detailed numerical analysis, we present the predictions 
for the inclusive and differential cross sections of the A°Z° associated production at the 
LHC. Sec. V contains a brief conclusion. For completeness, the relevant Feynman rules are 
collected in Appendix A, and the lengthy analytic expressions of the result of our calculation 
are summarized in Appendices B and C. 



II. LEADING ORDER CALCULATIONS 

The related Feynman diagrams which contribute to the LO amplitude of the partonic 
process b(pi)b(p 2 ) — * Z°(p 3 )A°(p4) are shown in Fig. 1. The LO amplitude in n = 4 — 2e 
dimensions is 



with 



M,S s) = 2im b ( AlF \ + A2F \ )v(p 2 )u(p 1 )p A .e p (p 3 ) , 
V s - m^o s - mjjo ) 

M§> = *-^v( P 2h 5 (A- A) tfpsXCV + C Al ")u{ Pl ) , 

M ( U) = hM^j An){Cv + CAl r> ){A _ A)7 5 m0oi)) 
u 

where 5^ is the color tensor (ii,?2 are color indices for the initial state quarks), /i r is a 
mass parameter introduced to keep the couplings dimensionless, s, t and u are Mandelstam 
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variables, which are defined as 

s = (pi + P2) 2 , * = (pi - P3) 2 , u = (p2 - Ps) 2 , 

A=i,2,3, ^i=i,2 and CV ;J 4 denote the coefficients appearing in the relevant h°(H°, A°)bb, 
h (H°)Z°A and Z°bb couplings, respectively, and their explicit expressions are given in 
Appendix A. 

In order to simplify the expressions, we further introduce the following Mandelstam vari- 
ables: 



t' — t — m| , u = u — m 2 z o ■ 



(1) 



After the n-dimensional phase space integration, the LO partonic differential cross sec- 
tions are given by 



2 ±.B 



d 2 a 



71 S f 



dt'du' 



t'u' — sm 



z° 



s 2 T(l -e) 
x5(s + t + u — m| - m 2 A o)^2\M 



Q(t'u — sm 2 z o)Q[s — (m z o + m A o) 2 ] 



B\2 



(2) 



with 



4m 2 A oS 



(s - m| - m 2 A o) 2 s 



m 



z° 



MF\ A 2 F 2 
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h" 



s — m 



+ A 2 3 (\C V \ 2 + \C A \ 2 ) 



+A 2 3 (\C V \ 2 + \C A \ 2 ) 
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2 2m| (l — e)(tu — m 2 z0 m A o) + su 2 



turn 
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m z0 u z 



(s-"ifco) (s-m| ) 



+2^(|CV| 2 - |C A | ; 



2(1 — e)m z o(tu — m 2 z0 m A o) + s(2m z0 m A0 — t-u) 
m 2 z0 tu 



,(3) 



where = (47r)~ 2+e and the O function is the Heaviside step function. J2\M B \ 2 is the LO 
squared matrix element of b(xipi)b(x 2 p 2 ) — > Z (p 3 )A°(p4), in which the colors and spins of 
the outgoing particles have been summed, and the colors and spins of the incoming ones 
have been averaged over. 

The LO total cross section at the LHC is obtained by convoluting the partonic cross 
section with the parton distribution functions (PDFs) G bb / p in the proton: 



a B = J dx 1 dx 2 [G b / p (x 1 , fif)G- b / p (x 2 , fif) + (xi <-> x 2 )]a B , 



(4) 



where /if is the factorization scale and a B is the Born level constituent cross section of 
b[x\ , p\)b{x2 r p < i) Z (p3)A (p^). Obviously, the above LO results in the DREG scheme are 
equal to the ones in the DRED scheme since the LO calculations are finite and free of any 
singularity. 



III. NEXT-TO-LEADING ORDER CALCULATIONS 

The NLO contributions to the associated production of A and Z° can be separated into 
the virtual corrections arising from loop diagrams of colored particles and the real corrections 
arising from the radiation of a real gluon or a massless (anti) quark. For both the virtual and 
real corrections, we will first present the results in the DREG scheme, and then compare 
them with the ones obtained in the DRED scheme. 



A. Virtual corrections 

The virtual corrections to bb — > A°Z° arise from the Feynman diagrams shown in Fig. 2 
and Fig. 3. They consist of self-energy, vertex and box diagrams, which represent the SM 
QCD corrections, arising from quarks and gluons, and supersymmetric QCD corrections, 
arising from squarks and gluinos. We carried out the calculation in t'Hooft-Feynman gauge 
and used the dimensional regularization in n = 4 — 2e dimensions to regularize the ultraviolet 
(UV), soft and collinear divergences in the virtual loop corrections. In order to remove the 
UV divergences, we renormalize the bottom quark masses in the Yukawa couplings and the 
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wave function of bottom quark, adopting the on-shell renormalization scheme 

Denoting m bQ and ip b o as the bare bottom quark mass and the bare wave function, re- 
spectively, the relevant renormalization constants Srrib, SZ b L and 5Z bR are then defined as 

m b0 = m b + 5m b , (5) 
ij b0 = (1 + 8Z bL )^ bL + (1 + 5Z bR )^ bR . (6) 

After calculating the self-energy diagrams in Fig. 2, we obtain the explicit expressions of all 
the renormalization constants as follows: 



5m b 
m b 



-g<^{3£oK,0,m 6 2 )-2 
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1=1 
a 



Tfl 

£ £i-^sin20 s (-l)\B c 



m b 



(ml m], mf.) f , 



5Z 6L = -^<7 F E(^i) 2 (^o + Bi) (0, mf. , m|), 



5Z, 



bR 



a s 
2^ 



i=l 
2 



CFE(^ 2 ) 2 (5o + 5i)(0,m?,m?), 



i=l 



where CV = f, B , B 1 are the scalar two-point integrals 26], are the sbottom masses, 



ffl» is the gluino mass, and R b is a 2 x 2 matrix shown as below, which is defined to transform 



the sbottom current eigenstates to the mass eigenstates |27]: 



R b 



R b 



>R 



cos 6~ b sin 6i 
- sin Q~ h cos Q~ h 



(7) 



with < &i < 7r, by convention. Correspondingly, the mass eigenvalues and (with 
nT>b < m b 2 ) are gi ven by 



1 ml 

bi 

m~ 
62 



R b Ml{R b )\ 



t 2 

ya b m b m\ R 



(8) 



with 



m 
in 



2 

a& 



M Q + m 6 + m l cos 2 /^( J 3L - e & sin2 8w), 
A b — fi tan /3. 



+ m 2 + m| cos 2/3eb sin 2 6*^, 



(9) 
(10) 

(11) 



Here, M~ is the sbottom mass matrix. Mq ^ and A b are soft SUSY breaking parameters 
and /j, is the higgsino mass parameter . /| L and e b are the third component of the weak 
isospin (i.e. —1/2) and the electric charge of the bottom quark b (i.e. —1/3), respectively. 
The renormalized virtual amplitudes can be written as 



M v = M v 



+ M C 



(12) 



Here, M unren contains the radiative corrections from the one-loop self-energy, vertex and 
box diagrams, as shown in Fig. 2, and M con is the corresponding counterterm. Moreover, 
M unren can be separated into two parts: 

g d 



M unren _ M « + M box(/3)^ 

ct=a (3= a 



(13) 



where a and (5 denote the corresponding diagram indexes in Fig. 2 and Fig. 3, respectively. 
They can be further expressed as 
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i=i 
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M box(f3) = Y j f b l OX(f)) M h 
1=1 
12 

j^unren "y y junren 

1=1 



(f? 



(14) 
(15) 
(16) 



where ff and fi° x ^ are the form factors, which are given explicitly in Appendix B, and the 
Mi are the standard matrix elements defined as 

Mi (2) = v(p 2 )u(p 1 )p 1{ 2).e(p 3 ), 
M 3 (4) = v{P2)PRu{pi)p 1 (2).e(p 3 ), 
M 5 = v(p 2 ) i> 3 ^(p 3 )u(pi), 
M 6 = v(p 2 ) i> 3 ^(p 3 )P R u(p 1 ), 
M 7 (8) = v(p 2 ) 3 PRu(pi)p m .e(p 3 ), 

Mg = V(p 2 ) /(Ps)m(Pi), 

Mio = v(p 2 ) tf{p 3 )P R u(pi), 

M 11(12) = v(p 2 ) ^ M (Pi)Pi(2)-e(P3)- (17) 

The counterterm contribution M con is separated into M con ( s \ M con( ^ and M con( - u \ i.e. the 
counterterms for s, t and u channels, respectively, which are given by 



j^con(s) _|_ j^con(t) _|_ j^con(u) 



M c 



s — m 



h° 



s — m 



H° 



5m b + ^-(5Z bL + 5Z 6i? ) 



(M x + M 2 ) , 



j^con(u) 



-iA 3 
t 

iA 3 
u 



m b 



5m b + —(5Z bL + 5Z bR 



[2(CV - C A )M l - AC V M 3 
-(C V -C A )M 5 + 2C V M 6 ], 



5m b + ^-(5Z bL + 5Z bR ) 



[2(C y + C A )M 2 - AC V M 4 



+{C V + C A )M 5 - 2C V M 6 ] . 



The 0{a s ) virtual corrections to the differential cross section can be expressed as 

1 Q(t u - sm 2 z o)Q[s - (m z o + m A i>y] 



dt'du' s 2 r(l-e) 



xS(s + t + u — 777/40 — rn%o) 2 Re 



Y J {M V M B *) , (18) 



where the renormalized amplitude M v is UV finite, but it still contains the infrared (IR) 
divergences: 

where 

A\ = -C F , A\ = -ho F . (20) 

Here, the infrared divergences include the soft divergences and the collinear divergences. 
The soft divergences are cancelled after adding the real emission corrections, and the re- 

rn 

maining collinear divergences can be absorbed into the redefinition of PDF 28], which will 
be discussed in the following subsections. Note that the coefficients A\ and A\ of the in- 
frared divergence terms are constants, similar to the Drell-Yan type processes. Needless to 
say that the SUSY QCD corrections do not generate infrared divergences, for squarks and 
gluinos are massive particles. 

In the above calculation, we have adopted the naive 75 prescription in the DREG scheme 
to calculate the A°Z° associated production rate. To cross check the above calculation, 
we shall also adopt the DRED scheme to carefully treat the 75 factor in the amplitude 
calculation. We shall show that the total inclusive rate is independent of the regularization 
scheme, though the individual contributions, from either virtual or real emission corrections, 
can be scheme-dependent. 

In the DRED scheme, 5Z b i and SZ b R remain unchanged, however, 5m b is different, and 

fSm b \ fSm b \ a s 

Up. [21) 



\m b J DREG \m b J DRED 4:71 

Similarly, the form factors are found to be different, and 

fr en dreg ~fr en DRED = ~^C F , for i = 1, 2...6, (22) 

DREG - DRED = ^ ^ • = 7, 8...12. (23) 

Thus, it is easy to obtain the following relations from the above results: 

M V D REG ~ M V D RED = ~^C F M B , (24) 

° V DREG - CT V DRED = ~Tr-C F a B + 0(0%). (25) 

Z7T 
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B. Real gluon emission 



The Feynman diagrams for the real gluon emission process b(pi)b(p2) — > Z°(ps)A°(p4) + 
g(p§) are shown in Fig. 4. 

The phase space integration for the real gluon emission will produce infrared singularities, 
which can be either soft or collinear and can be conveniently isolated by slicing the phase 
space into different regions defined by suitable cut-offs. In this paper, we use the two- 
cutoff phase space slicing method which introduces two small cut-offs to decompose the 
three-body phase space into three regions. 

First, the phase space can be separated into two regions by an arbitrary small soft cut-off 
8 a , according to whether the energy (E 5 ) of the emitted gluon is soft, i.e. E 5 < 8 s yfs/2, or 
hard, i.e. E 5 > 5 s y/s/2. Correspondingly, the partonic real cross section can be written as 

a R = c s + * H , (26) 

where a s and a H are the contributions from the soft and hard regions, respectively. a s 
contains all the soft divergences, which can be explicitly obtained after analytically in- 
tegrating over the phase space of the emitted soft gluon. Second, in order to isolate 
the remaining collinear divergences from a H , we should introduce another arbitrary small 
cut-off, called collinear cut-off 8 C , to further split the hard gluon phase space into two 
regions, according to whether the Mandelstam variables satisfy the collinear condition 
— 5 c s < ui t 2 = (pi,2 — P5) 2 < or not. Thus, we have 

H ~HC 1 ±MC 



a = a 



+ a HC , (27) 



where the hard collinear part a HC contains the collinear divergences, which can be explicitly 
obtained after analytically integrating over the phase space of the emitted collinear gluon. 
The hard non-collinear part a HC is finite and can be numerically computed using standard 
Monte-Carlo integration techniques [3^, and can be written in the form: 

da WU =-Y / \M bh \ 2 dT 3 . (28) 

Here, dT^ is the hard non-collinear region of the three-body phase space. 

In the next two subsections, we will discuss in detail the soft and hard collinear gluon 
emission. 
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1. Soft gluon emission 



In the soft limit, i.e. when the energy of the emitted gluon is small, with E$ < 5 s ^/2, 
the matrix element squared J2\M R \ 2 for the process b(pi)b(p 2 ) —> Z (p 3 )A (p^)g(p 5 ) can be 
simply factorized into the Born matrix element squared times an eikonal factor $ e jfe: 



J2\M R (bb - A°Z° + g)\ 2 ^ (4vra^)E|M B r$ e , fc , 



(29) 



where the eikonal factor <& ei k is given by 



(Pi 'Ps)(P2'Ps)' 

Moreover, the phase space in the soft limit can also be factorized as 



(30) 



dT 3 (bb -> A U Z U + g) s ^ dT 2 {bb -> A"Z»)dS. 
where dS is the integration over the phase space of the soft gluon, which is given by [2 

rSsy/a/1 



(31) 



dS 



1 fSsVs/2 

/ dE 5 El- 2 ^n 2 _ 2e . 

it r Jo 



2(27r)--« Jo 

Hence, the parton level cross section in the soft region can be expressed as 



(32) 



a s = (47ra s/ ^) / dT 2 J2\M B \ 2 I dS$ eik . 



(33) 



Using the approach of Ref. [29], after analytically integrating over the soft gluon phase space, 
Eq. (j3^|) becomes 



a = a 



a s T(l — e) / Att^ 



27rr(l-2e)V s 



(34) 



with 



A s 2 = 2C F , A{ = -4C F In 8 S , A S Q = AC F In 2 S s . 



(35) 



2. Hard collinear gluon emission 

In the hard collinear region, i.e. E 5 > 5 s y/s/2 and — 5 c s < ui i2 < 0, the emitted hard 
gluon is collinear to one of the incoming partons. As a consequence of the factorization 



theorems 



the squared matrix element for bb — > A°Z° + g can be factorized into the 
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product of the Born squared matrix element and the Altarelli-Parisi splitting function for 
6(6) -> b(b)g HQ, i-e. 



Y}M R (bb^ A°Z° + g)f 



+ 



(36) 



1 2 collinear 

( |-n s //; i > | w - | i 

ZU\ ZU2 

where z denotes the fraction of incoming parton 6(6) 's momentum carried by parton 6(6) 
with the emitted gluon taking a fraction (1 — z), and Py (z, e) are the unregulated splitting 
functions in n = 4 — 2e dimensions for < z < 1, which can be related to the usual 



(37) 
(38) 



Altarelli-Parisi splitting kernels |22j as Pij(z,e) = Pij(z) + tP[Az). Explicitly 



P hh (z) = P- b i{z 



Pl b (z) = P&z) 



1 4- t 2 3 
C Fj ^— + C F -5(l-z), 

-C F (l-z)+C F ^5(l-z). 



Moreover, the three-body phase space can also be factorized in the collinear limit, and, for 
example, in the limit — 5 c s < u\ < it has the following form 

(4tt) 



dT 3 (bb -> A°Z° + g) 



dT 2 (bb -> A°Z"; s' = zs) 



167r 2 rfl - e 



-dzdu^z - l)wi]- e . (39) 



Here, the two-body phase space should be evaluated at the squared parton-partpn energy 
zs. Thus, the three-body cross section in the hard collinear region is given by 



da HC = a B 



a s r(l-e) ,4tt^, 



L27rr(l-2e) v s 
+P55(^,e)Gb /p (xi/2;)G';,/p(x2) + (x x <-> x 2 ) 



P bb (z,e)G b /p(x 1 /z)G- b/p (x 2 ) 
dz , 1 — z , 



e dx\dx2 ) 



(40) 



where G b (p\/ p (x) is the bare PDF. 



C. Massless (anti) quark emission 

In addition to the real gluon emission, a second set of real emission corrections to the 
inclusive production rate of pp — > A°Z° at the NLO involves the processes with an additional 
massless (anti) quark in the final states: 

gb -> bA°Z°, gb -> L4°Z° . 

The relevant Feynman diagrams for massless (anti) quark emission (the diagrams for the 
antiquark emission are similar and omitted here) are shown in Fig. 5 
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Since the contributions from the real massless (anti) quark emission contain the initial 
state collinear singularities, we also need to use the two-cutoff phase space slicing method 



29| to isolate those collinear divergences. Because there is no soft divergence in the splitting 
of g — > bb, we only need to separate the phase space into two regions: the collinear region 
and the hard non-collinear region. Thus, according to the approach shown in Ref. |2^ |. the 
cross sections for the processes with an additional massless (anti) quark in the final states 
can be expressed as 



da 



add 



a c (a(3 — > A°Z° + X)[G a / p (x 1 )G l 3/ p (x 2 ) + (xi <-> x 2 )]dxidx 2 

(a=g,/3=b,b) 



+& 1 



a s r(l - e) 47r^, 



27rr(l-2e) v s 



+G b/p (x 1 )P- bg (z,e)Gg /p (x 2 /z) + (x 1 «-> x 2 ) 



Pbg( z , t)Gg/ p (Xi/ z)G b / p (x 2 ) 

dz (\ — z\~ £ 



z \ z 



dx±dx 2 , 



(41) 



where 



P bg (z)=P- hn (z) = -\z 2 + (l-z) 



hg 



n g (z)=Pi g (z) = -z(l-z). (42) 



The first term in Eq. (}4~T|) represents the non-collinear cross sections for the two processes, 
which can be written in the form: 



da c = ^£\M^\ 2 dY„ 



(43) 



where a and (3 denote the incoming partons in the partonic processes, and dT 3 is the three 
body phase space in the non-collinear region. The second term in Eq. (|41jl represents the 
collinear singular cross sections. 



D. Mass factorization 



As mentioned above, after adding the renormalized virtual corrections and the real cor- 
rections, the partonic cross sections still contain the collinear divergences, which can be 
absorbed into the redefinition of the PDF at NLO, in general called mass factorization 
This procedure in practice means that first we convolute the partonic cross section with the 
bare PDF G a / P (x), and then rewrite G a / P (x) in terms of the renormalized PDF G a / P (x, /^/)in 
the numerical calculations. In the MS scheme, the scale dependent PDF G a / P (x, ///) is given 
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by 



G, 



G 



1 dz 

-P a p{z)G p/p {x/z). (44) 

x Z 



After replacing the bare PDF by the renormalized MS PDF and integrating out the collinear 
region of the phase space defined in the two-cutoff phase space slicing method [29], the 



resulting sum of Eq. (|471|) and the col 
remaining O collinear contribution as 



inear part (the second term) of Eq. (|4Tj) yields the 



29j: 



a 



coll 



a 



O'.s 



r(l-e) 



+ E 



27rr(l-2e)V s 
Af(a — > ag) 



{G b /p(xi, fif)G b / p (x 2 , Hf) + G b / P (xi, fif)G b / p (x 2 , /i/) 



+ A s c (a -> ag) 



ce=b,b 

+ (xi <-> x 2 )}dxidx 2 , 



G b / P (xi, fif)G b / p (x 2 , /i/) 



(45) 



where 



A[ c (b -> bg) = A[ c (b -> bg) = C F {2lnS s + 3/2), 



A S C = ATH-), 

x,fJif)= I ^ —Gf3/ P (x/y,fXf)P al3 (y) 



G 



a(=b,b)/pl J ') / 



P=g,a 



(46) 
(47) 

(48) 



with 



y 



(49) 



The NLO total cross section for pp — > A°Z° in the MS factorization scheme is obtained 
by summing up the Born, virtual, soft, collinear and hard non-collinear contributions. In 
terms of the above notations, we have 



a 



NLO 



dx\dx 2 { 



G b / P (x 1} /j,f)G b/p (x 2 , jj f ) + (xi <-> x 2 ) 



v cr +cr +cr + a J) + G 



coll 



+ X/ / dx\dx 2 

{a=g,/3=b,b) 



G a /p(xi, flf)G/3/ p (x2, ///) + (Xi <-> X 2 ) 



a c {a(3^ A U Z° + X).(50) 



o 



We note that the above expression contains no singularities, for 2A\ + A\ = and 2A] 7 + 
+Af (b -> 6#)+Af (6 -> 6flf) = 0. Namely, all the 1/e 2 and 1/e terms cancel in a NL °. The 
apparent logarithmic S s and 5 C dependent terms also cancel with the the hard non-collinear 
cross section a HC after numerically integrating over its relevant phase space volume. 
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E. Real emission corrections and NLO total cross sections in the DRED scheme 



In the end of Sec. Ill A, cf. Eqs. (|2T|l - (j25j) . we discussed the results of virtual corrections 
in the DRED scheme. Here, we examine the real emission corrections and the NLO total 
cross section in the DRED scheme and compare them with those obtained in the DREG 
scheme. We find that the contributions from soft gluon emission remain the same, while the 
ones from hard collinear gluon emission and massless (anti) quark emission are different due 
to the difference in the parton splitting functions and the perturbative PDFs. 

First, the splitting functions in the DRED scheme contain no e parts, so that 

Pij( z , c)dred = Pij( z )- (51) 
Thus, from Eqs. (|45j) and 1)51)1. we find the difference 

(T CM DREG ~ CT CoU DRED = --^ / a B {Y, / — Gp/^/y, fl f )P' b p(y)G b/p (x 2 , flf) 

ATI J a Jxi y 

+ 22 — G/3/ p (x 2 /y, fif)P i /3 (y)G b / p (xi, fif) + (x x <-> x 2 )}dx 1 dx 2 . (52) 

p J %2 y 

Seco ndly , the perturbative PDFs defined in the DRED and DREG schemes are different, 
and |34j: 

Ga/p\X, 

^})DREG - G a / p (x, Hf) D RED = 7T / P'ap(%/ y)G a / p (x, ^f) DREG- (53) 

Z7T p Jx y 

After substituting them into the formula for calculating the Born level cross sections, cf. 
Eq. P)l. we find the difference arising from the perturbative PDFs, at the 0(a s ) level, as: 



& DREG — & DRED — 7T~ / ° i / / 

2tt J ^ J x x y 



{^2 I —Gf3/ p (xi/y, fif)DREDP'bf3(y)Gl/ p (x2, fif)DRED 
+ S / G p /p (x 2 /y, Hf)DREDP'lp(v)Gb/p{x\> l^f) DRED + {xi «-> X 2 )}dx 1 dx 2 . (54) 

p Jx2 y 

Except the upper limit of the integral over y, the two expressions in Eqs. 1)5*2*)) and 1)51)1 are 
the same. After substituting Eqs. ()52jl . ()5*4*)l and ()2*5)l into Eq. 1)50)1. we find the relation 
between the two NLO total cross sections, separately calculated in the DREG and DRED 
schemes, as follows: 

a NLO D REG ~ a NLO DRED = f a B {j: f -Gp/^/y, H )P' b p(y)G b/p (x 2 , H ) 
ri dy 

+ Y. —Gp/ P (x 2 /y,n f )P'- b p(y)G b/p (x 1 ,fXf) + [x x «-> x 2 )}dx 1 dx 2 

p Ji - s s s bi3 y 

-^C F a B DRED + 0{a% (55) 
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Using the explicit expressions of the e parts of the splitting functions P', cf. Eqs. and 
P2jh we find 



^DREG = a DRED + U { a s)- ( 5 °) 

As expected, both schemes yield the same NLO total cross sections, up to 0(a 2 ). 

F. Differential cross sections in transverse momentum and invariant mass 

In this subsection ,we present the differential cross section in the transverse momentum 
of Z° and A bosons, respectively, and the invariant mass of the A°Z° pair. Using the 



351 ]. the differential distribution of the transverse momentum [px) 



notations defined in Ref. 
and rapidity (y ) of Z° boson for the processes 

P(pi) +P(P2) - Z°(p 3 ) + A ( Pi )[+g(p 5 )/b(p 5 )/b(p 5 )] (57) 

is given by 

f ° A = %PtS I dx i I dx 2 x 1 G a / p (x 1 , Hf)x 2 Gp/ p (x 2 , ///) ^ f."^ , (58) 
dpxdy q /3 i 2 

where yS is the total center-of-mass energy of the collider, and 

2 T 2 U 2 2 1 , /^2x 

Pt = ~S mz °' V= 2 W 2 '' 

—T 2 — m 2 7 o+m 2 A o _ -xiU 2 - m| + m 2 A0 
x i = ^T77 ' x 2 = ^—7^ (59) 

with T 2 = (p 2 — P3) 2 — m| and C/ 2 = — P3) 2 — m| - The limits of integral over y and px 
are 

- y ma * (pt) <y< y max (pt) , < PT < V T X , (60) 



with 



y max (pT) = arccoshf 



S + m| - m 2 A o 



^2JS(p 2 T + m 2 z0 ] 



1 



pT x = ^^(S + mlo - m\ o y - Am 2 z0 S . (61) 
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The differential distribution with respect to pt and y of A is similar to the one of Z°. The 
differential distribution with respect to the invariant mass Maz is given by 

da 2M/ 



1 az dCjiz 



dM AZ 



E 



S a dr 

a,p 



where dC , AZ is the parton luminosity, defined as: 



OL0 



J pa/3 

dr 



1 dx 

x 



Ga/p(x,fJ.f)G 0/p (T/x,fif) 



with 



(62) 



(63) 



M AZ = \J(E 3 + E 4 ) 2 - (p| + pt) 2 > (m A o + m z o), 
r = M\ z /S. 



(64) 
(65) 



IV. NUMERICAL RESULTS 



In the numerical calculations, we used the following set of SM parameters 

a ew (m w ) = 1/128, m w = 80.419 GeV, m z = 91.1882 GeV . 
m t = 178 GeV, a s (M z ) = 0.118. 



(66) 



The running QCD coupling a s (Q) is evaluated at the two-loop order |37j, and the CTEQ6M 
PDFs j3| is used throughout this paper to calculate various cross sections, either at the LO 
or NLO. As to the Yukawa coupling of the bottom quark, we shall first use the MS bottom 
quark mass, m&(m&) = 4.25 GeV, to evaluate the event rate, then compare it with the one 
calculated using the QCD improved running mass to reduce the higher order QCD radiative 
corrections, therefore improve the perturbative calculations. The QCD improved running 
mass rrib{Q), evaluated by the NLO formula 3^, is: 



m b (Q) = U 6 (Q,m t )U 5 (m t ,m b )m b (m b ) , 



(67) 



where the evolution factor Uf is 
U f {Q 2 ,Qt) 



oi 8 {Q\ 



12 



33-2/' 



jif) 



1 a s (Qi) - a s (Q2) j( f ) 
8982 - 504/ + 40/ 2 



3(33 -2ff 



(68) 
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and / is the number of the active light quarks. For comparison, we list the QCD improved 
running bottom quark mass in Table |l] for various energy scale Q. 

For large tan /3, the SUSY threshold correction to the bottom quark Yukawa couplings 
could be large, and it can be resummed by making the following replacement in the tree-level 
couplings to improve the perturbation calculations j^]: 

(n\ rn b (Q) /an s 

mbiQ) ~* l + A mb (Q = M SU s Y y m 



2a s (Q = M SUSY ) ^ ^ QT( , h 



Am b = — MgH tan (5 1 (m~ bi ,m~ b2 , Mg) + j^/iAttan/^m^, m i2 , fi) 



9 Z 



16vr 2 

1=1 



fiM 2 tan/3 ^ (Ri) 2 I(m k ,M 2 ,v) + -(Ri) 2 I(m k ,M 2 ,ti) , (70) 



where 



h = , (72) 

y2mvi/ sin p 

and i?* and are the rotation matrices for defining the mass eigenstates of t and b, respec- 
tively. We set Msusy in Am& to m g in our numerical calculations. Needless to say that 
when using the running bottom quark Yukawa coupling to evaluate cross sections, we shall 
subtract the corresponding (SUSY-) QCD corrections at the order a s from the renormaliza- 
tion constant 8m b to avoid double counting in perturbative expansion of the strong coupling 
constant. 

The values of the MSSM parameters taken in our numerical calculations were constrained 
within the minimal supergravity scenario (mSUGRA) 0|, in which there are only five free 
input parameters at the grand unification (GUT) scale. They are mi , mo, Aq, tan/3 and 
the sign of fi, where mi, mo, A) are ; respectively, the universal gaugino mass, scalar mass 
and the trilinear soft breaking parameter in the superpotential. Given those parameters, 
all the MSSM parameters at the weak scale are determined in the mSUGRA scenario by 



using the program package SUSPECT 2.3 



421 ] . In particular, we used the running Higgs 



masses at the m z scale, defined in the modified dimensional reduction (DR) scheme, which 
have included the full one-loop corrections, as well as the two-loop corrections controlled 
by the strong gauge coupling and the Yukawa couplings of the third generation fermions 
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Q (GeV) 


250 


500 


750 


m b (Q) (GeV) 


2.68 


2.55 


2.49 



TABLE I: The QCD improved running bottom quark mass, evaluated at Q 
750 GeV. The MS bottom quark mass is taken to be mb(nib) = 4.25 GeV. 



250,500, and 



421 |43(. In our numerical calculations, we used the two- loop renormalization group equa- 
tions (RGEs) presented in that program for calculating all the gauge couplings, the (third 
generation) Yukawa couplings and the gaugino masses, while using one-loop RGE for the 
other supersymmetric parameters. In the following, we shall present our numerical studies 
based on the five sets of SUSY input parameters listed in Table ITU which are consistent with 
all the existing experiment data j^fj]. We will also vary tan/5, mo and Aq to examine their 
effects to various cross sections. For completeness, we also show the relevant SUSY output 
parameters in Table ILTT1 The QCD plus SUSY-QCD and SUSY-EW improved bottom quark 
running mass are listed in Tabl dlVl which should be compared with those given in Table |TJ 
in which only QCD running effect is included. For comparison, the QCD plus SUSY-QCD 
improved bottom quark running mass are separately listed in TablJVl 



set number 


m (GeV) 


mi (GeV) 


A, (GeV) 


tan (3 


sign(At) 


1 


150 


180 


300 


40 


+ 


2 


150 


400 


300 


40 


+ 


3 


200 


160 


100 


40 




4 


250 


160 


100 


40 




5 


400 


160 


100 


40 





TABLE II: Five sets of SUSY input parameters studied in this paper, within the mSUGRA 



scenario. 



As for the renormalization and factorization scales, we always chose [l t = m av 
m z o)/2 and [if = m av , unless specified otherwise. 



(m A o + 
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m l m ( GeV ) 

°1(2) 




m g {GeV) 


m A°(h°,H°)(GeV) 


A(b)(GeV) 


/i(GeV) 


a 


COS 97,7, 


1 


374.6(429.1) 


339.7(457.7) 


457.0 


223.8(107.5,223.9) 


-256.5(-275.8) 


235.3 


-0.032 


0.97(0.74) 


2 


764.3(822.0) 


673.6(833.8) 


940.0 


458.3(115.5,458.3) 


-607.5(-750.9) 


498.1 


-0.027 


0.47(0.71) 


3 


314.3(395.1) 


305.5(425.1) 


416.8 


133.7(106.7,134.1) 


-263.9(-303.9) 


-224.6 


-0.143 


0.62(-0.69) 


4 


330.8(408.6) 


317.0(434.3) 


419.9 


155.0(107.2,155.3) 


-262.7(-303.6) 


-228.8 


-0.086 


0.60(0.71) 


5 


396.1(467.5) 


363.9(476.1) 


431.9 


233.0(108.4,233.2) 


-261.0(-304.9) 


-249.4 


-0.043 


0.54(0.79) 



TABLE III: The SUSY output parameters used in our numerical calculations, corresponding to 
the five sets of SUSY input parameters listed in Table [H] 



set number 


1 


2 


3 


4 


5 


m b (Q = 250, Msusy = m g )(GeV) 


2.35 


2.41 


3.18 


3.16 


3.10 


m b (Q = 500, Msusy = m s )(GeV) 


2.24 


2.29 


3.03 


3.01 


2.96 


m b (Q = 750, Msusy = m g )(GeV) 


2.18 


2.23 


2.95 


2.93 


2.88 



TABLE IV: The QCD plus SUSY-QCD and SUSY-EW improved bottom quark running mass for 
the five sets of SUSY inputs listed in Table ITT1 evaluated at Q = 250, 500 and 750 GeV. 



set number 


1 


2 


3 


4 


5 


m b (Q = 250, M SUSY = m- g )(GeY) 


2.15 


2.14 


3.71 


3.66 


3.53 


m b (Q = 500, Msusy = m g )(Ge\) 


2.04 


2.04 


3.53 


3.49 


3.36 


m b (Q = 750, Msusy = m g )(GeV) 


1.99 


1.98 


3.44 


3.39 


3.27 



TABLE V: The QCD plus SUSY-QCD improved bottom quark running mass for the five sets of 
SUSY inputs listed in Table ILT1 evaluated at Msusy = rn§, and Q = 250, 500, and 750 GeV. 

A. LO total cross section 

In Fig. 6 and Fig. 7, we first compare the LO total cross sections of pp — > A°Z° via bb 
annihilation with the ones via gluon fusion and Drell-Yan processes, respectively. Here, we 
use the MS bottom quark mass m b (m b ) = 4.25 GeV, without including the effect from QCD 
running. Our numerical results are different from the ones presented in Ref. ^j], because 
the updated SUSY parameters are used instead of the earlier input parameters used in 
Ref. jl8| which have already been ruled out by recent experiments. As shown in Figs. 6 
and 7, the LO total cross sections via bb annihilation and Drell-Yan processes increase with 
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tan /3, while the ones via gluon fusion process are relatively larger for low and high values 
of tan/?, but become smaller for intermediate values of tan/3. Moreover, all the LO rates 
decrease when m^o increases. Figs. 6 and 7 also show that in most of the chosen parameter 
range, bb contributions are much larger than the ones from gluon fusion and Drell-Yan 
processes, especially for large tan/3 and small m^o, where the total cross sections from the 
bb contributions can reach a few hundred fb. 

B. Cutoff dependence 

In Fig. 8, we show the dependence of the NLO QCD predictions on the two arbitrary 
theoretical cutoff scales 8 S and 8 C , introduced in the two-cutoff phase space slicing method, 
where we have set 5 C = 8 s /50 to simplify the study and used QCD plus SUSY improved 
bottom quark Yukawa coupling. The NLO total cross section can be separated into two 
classes of contributions. One is the 2 — > 2 rate contributed by the Born level, and the 0(a s ) 
virtual, soft and hard collinear real emission corrections, denoted as a B , a , cr s and o co11 
in Eq. ()50|). Another is the 2 — + 3 rate contributed by the 0(a s ) hard non-collinear real 
emission corrections, denoted as a HC and a c in Eq. (jSUJ). As noted in the previous section, 
the 2^2 and 2 — > 3 rates depend individually on 5 S and 5 C , but their sum should not 
depend on any of the theoretical cutoff scales. This is clearly illustrated in Fig. 8 for two 
different sets of SUSY parameters. We find that o^lo is almost unchanged for 5 S between 
5 x 10~ 5 and 10~ 2 , which is about 200 fb and 28 fb, respectively for the two different sets of 
SUSY parameters. Therefore, we take 5 S = 10 -4 and 5 C = 5 s /50 in the numerical calculations 
below. 

C. m A o dependence 

In Fig. 9, we show the total cross sections of pp — > A°Z° at the LHC as a function of m^o 
for tan/3 = 10 and 40, respectively, assuming mi = 160 GeV, A = 100 GeV, and /i < 0. 
We considered the LO total cross sections in three different cases, i.e. using (I) MS bottom 
quark mass at the scale m^, (II) QCD improved bottom quark running mass at the scale 
m^o, and (III) QCD plus SUSY improved bottom quark running mass at the scale m^o, 
respectively. We also considered the NLO total cross sections for the cases of (II) and (III). 
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Fig. 9 shows that the LO and NLO total cross sections get smaller with the increasing m A o, 
and the results for tan (3 = 10 in Fig. 9(2) are much smaller than the ones for tan/3 = 40 
in Fig. 9(1). For small m A o (< 160 GeV) the LO total cross sections in Fig. 9(1) can be 
larger than 100 fb. The contributions from the QCD running m& mass effects and the SUSY 
improved rrif, corrections are significant, for example, in Fig. 9(1) when m A o ~ 155 GeV and 
tan P = 40, the LO total cross sections are about 270 fb, 120 fb and 185 fb for the three cases, 
respectively. Moreover, Fig. 9 shows that the NLO QCD corrections can either enhance or 
suppress the total rate, and the 0(a s ) contribution is in general a few tens percent of the 
total rate, as described below. Define the K-factor as the ratio of the NLO to LO total cross 
sections, calculated using the CTEQ6M PDFs. We shown in Fig. 10 the dependence of the 
K factor on m A o for A°Z° production, based on the results of case (II) in Fig. 9. Namely, the 
QCD improved bottom quark running mass is used for calculating the total cross section at 
the LO and the NLO. Fig. 10 shows that in general the K factor becomes smaller with the 
increasing m A o. For example, the curve (a) in Fig. 10(1) shows that when m A o varies from 
108 GeV to 900 GeV, the K factor varies from 1.72 to 0.82, and the curve (a) in Fig. 10(2) 
shows that when m A o varies from 235 GeV to 860 GeV, the K factor varies from 0.91 to 0.68. 
The contributions to the K factors, shown as curve (a), in both Fig. 10(1) and Fig. 10(2) 
come from the pure QCD corrections, shown as curve (b), and SUSY QCD corrections, 
shown as curve (c). The former includes both the virtual and real emission contributions 
originated from pure QCD corrections, while the latter consists of only virtual corrections. 
As expected, the K-factor contributed by the pure QCD corrections is under controlled, 
of a few tens percent, when the QCD improved bottom quark running mass is used to 
evaluate the Yukawa coupling of the bottom quark. On the other hand, the SUSY QCD 
corrections could become large as m A o decreases, especially for large tan (3. For example, in 
Fig. 10(1), for tan/5 = 40, when m A o ~ 108 GeV, the K factor of SUSY QCD corrections 
is about 0.8 which dominates the overall K factor. Hence, to improve the convergence of 
the perturbation calculations in the case of large tan f3, we could use the SUSY improved 
bottom quark running mass to evaluate the Yukawa coupling of bottom quark. More on 
SUSY QCD corrections will be discussed below. We have also examined the contributions 
from the box diagrams shown in Fig. 3. The pure QCD box diagram contribution, arising 
from Fig. 3(a) and Fig. 3(c), is ultraviolet finite but not infrared finite. For tan/5 = 40 the 
finite part of the pure QCD box diagram contribution becomes more important for large 



21 



m A o, and its effect is to decreases the total rate. On the contrary, the SUSY QCD box 
diagram contribution, arising from Fig. 3(b) and Fig. 3(d), is free of any singularity, and is 
small numerically. 

Fig. 11 shows the dependence of the K factors on m A o for A°Z° production, based on 
the results of case (III) in Fig. 9. Namely, the QCD plus SUSY improved bottom quark 
running Yukawa coupling is used for calculating the total cross section at the LO and the 
NLO. Generally, the K factor decreases with m A o. For example, for tan/3 = 40, when m A o 
varies from 108 GeV to 900 GeV, the K factor corresponding to curve (a) ranges from 0.98 
to 0.61, which contains two parts: the pure QCD corrections, shown as curve (b), and SUSY 
QCD corrections, shown as curve (c). As compared to the results in Fig. 10(1), we find that 
the SUSY QCD correction, shown as curve (c), has been largely suppressed. For instance, 
the K factor of SUSY QCD corrections drops from 0.8, in Fig. 10(1), to 0.05, in Fig. 11, for 
m A o ~ 108 GeV, while the other SUSY parameters are identical in both calculations. This 
is because using the SUSY improved running m h to evaluate the LO cross section, we have 
already included the dominant NLO SUSY QCD corrections. Therefore, we shall use the 
QCD plus SUSY improved bottom quark running mass in the following numerical analysis 
for both the LO and NLO calculations, unless specified otherwise. 

D. SUSY QCD corrections in heavy mass limit 

It is instructive to examine the results of Figs. 10(1) and Fig. 11 in the heavy mass limit, 
where all the SUSY mass parameters except m A o are of the same size and tend to be heavy, 
i.e. Mq = M b = n = A b = M- g = M SU sy > m Z - In the heavy mass limit, the SUSY QCD 
box diagram contribution, arising from Fig. 3(b) and Fig. 3(d), is suppressed by powers of 
Msusy and can be neglected. This is confirmed by our numerical calculation which shows 
that the SUSY QCD box contribution is generally below 0.1% of the total rate. Hence, we 
shall examine the effect of SUSY QCD corrections in the heavy mass limit to the virtual 
diagrams shown in Figs. 2(a)-(g), and compare the analytical result with our numerical 
calculations. 

Since our aim is to examine the NLO SUSY QCD effect in this part of study, we shall 
use the LO bottom quark Yukawa coupling (with m& = 4.25 GeV) to evaluate the relevant 
tree level vertices. Keeping only terms at 0(a s ) that are not suppressed by negative powers 
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of heavy mass Msusy in the heavy mass limit, the one loop SUSY QCD correction to the 
individual diagram in Fig. 2 yields the following corrections. After stripping off the Born 
level matrix element (including all the vertex and propagator factors), the multiplicative 
factor of the s-channel diagram with the h° propagator, cf. Fig. 2(a), is given by 

F Hh o = -^(l + cota) (73) 



where a is the mixing angle of the two CP-even Higgs bosons |2j. Note that Eq. (J73*)l is 
in agreement with the one shown in Ref. j^J. Similarly, the multiplicative factor of the 
s-channel diagram with the H° propagator, cf. Fig. 2(a), is given by 

F (a) „o = -^(l-tana). (74) 

The multiplicative factor of the t and u-channel diagrams, cf. Fig. 2(c) or (d), is given by 

F (c)A o = F (d)A o = -yf^(l + cot P) . (75) 

The multiplicative factor for the sum of Figs. 2(b) and 2(g) is zero. This is because after 
adding the wavefunction renormalization factor for the external bottom quark line, the 
renormalized Zbb vertex vanishes in the heavy mass limit. (Again, we have dropped any 
term that is suppressed by negative powers of the heavy mass scale Msusy-) Similarly, the 
multiplicative factor for the sum of Figs. 2(e) and 2(f) is zero. 

Given the above multiplicative factors, we can calculate the SUSY QCD correction to the 
total cross section for pp — > A°Z° production at the LHC, and compare it with the complete 
numerical calculation described in Sec. III. The results in the heavy mass limit are shown in 
Fig. 12, which show that the agreement becomes better for larger value of Msusy- Hence, 
this provides a consistent check on our complete numerical calculations. 



E. tan/3 dependence 

In Figs. 13 (1) and 13 (2), the total cross sections for pp — > A°Z° at the LHC are 
plotted as a function of tan/? for two representative values of mi and tuq, respectively. In 
Fig. 13(2), when tan/3 ranges between 4 and 40, m A o varies from 330 GeV to 223 GeV, and 
from 660 GeV to 458 GeV for mi = 180 GeV and 400 GeV, respectively. From Fig. 13(2) we 
can clearly see that the LO and NLO total cross sections are enhanced with the increasing 
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tan/3 and decreased with the increasing mi. For large tan/3 (> 40) and mi = 180 GeV, the 
LO and NLO total cross sections can be over 30 fb. The features in Fig. 13(1) are similar 
to the ones in Fig. 13(2), but in general the total cross sections are larger than later. For 
example, for large tan/?(> 40) and mo = 200 GeV, both of the LO and NLO total cross 
sections can reach about hundreds of fb. 

Fig. 14 shows the dependence of the K factors on tan/?, based on the results in Fig. 13, 
where the K factor increases with the increasing tan/3. For the results of Fig. 13(1), the 
K factor varies from 0.69 to 0.90 and from 0.65 to 0.92 for m = 200 GeV and 400 GeV, 
respectively. For the results of Fig. 13(2), the K factor varies from 0.70 to 0.74 and from 
0.62 to 0.63 for mi = 180 GeV and 400 GeV, respectively. 

F. fJ-r/^f dependence 

Fig. ITB1 shows the dependence of the total cross sections for pp — > A°Z° production at the 
LHC on the renormalization scale (// r ) and the factorization scale (///), with \i r = \x$. The 
case (1) is for /x < 0, and the case (2) is for \i > 0. In both cases, the scale dependence of 
the NLO total cross section is smaller than that of the LO cross section. For example, the 
LO cross sections vary from 65 fb to 261 fb and 25 fb to 45 fb when /j r = /zj ranges between 
0.1 m av and 10m av , while the NLO ones vary from 230 fb to 232 fb and 38 fb to 39 fb, in the 
case (1) and (2), respectively. Here, the QCD plus SUSY improved bottom quark Yukawa 
coupling is used. For comparison, we also show the results of other two calculations. The 
case (3) is similar to the case (1), but in (3) the pure QCD running bottom quark mass 
is used instead. The case (4) is similar to the case (1), but in (4) the contribution from 
the SUSY-EW correction in the running bottom quark Yukawa coupling is not included, 
namely, only the pure QCD and SUYSY-QCD corrections are included. 

To further investigate the scale dependence in case (1), with /i < 0, we study the scale 
dependence of the total cross section on the renormalization scale (fi r ) and the factorization 
scale (/i/) seperatedly in Fig.UHl Here, the QCD plus SUSY improved bottom quark Yukawa 
coupling is used. We find that in either case, whether we fixed \i r and let /i/ vary, or vice 
versa, the NLO rate is less dependent on the scale than the LO rate. 

Hence, when applying the usual prescription to estimate the scale dependence, i.e. varying 
the scale around m av by a factor of 2, the NLO cross sections vary by around 10% to 20%, 
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cf. Figs. M and DH 



G. PDF uncertainty 

To estimate the uncertainties in the total cross sections due to the uncertainty of PDFs, 
we take the 41 sets of CTEQ6.1 PDFs to calculate the LO and NLO rates [46]. As shown 
in Fig. El the LO result of using the CTEQ6M PDF lies between the maximum (a max ) 
and minimum (<7 m j n ) LO rates. The NLO total cross sections are then calculated using 
three different PDF sets, one of which is CTEQ6M, the other two are the ones that give 
the maximum and minimum LO rates, respectively. The total cross sections for pp — > A°Z° 
production at the LHC, as a function of the trilinear coupling Aq, for the above mentioned 
PDFs are shown in Fig. El where we have used the QCD running mass to evaluate the 
bottom quark Yukawa coupling. It turns out that the PDF uncertainties (defined here as 
i(cmaa: — &min)) / (&max + °"mm)) in the LO and NLO total cross sections are about the 
same, when the QCD running m& is used. For example, when Aq = 100 GeV, the PDF 
uncertainties are ±2.9% at the LO, and ±3.0% at the NLO, respectively. 

Fig. ITS1 shows the PDF uncertainties (defined here as the Eq. (3) in Ref. 0|) in the LO 
and NLO total cross sections for pp — > A°Z° production at the LHC, as a function of m^o. 
Here, we also used the QCD running mass to evaluate the bottom quark Yukawa coupling. 
It turns out that the NLO rate has a slightly larger uncertainty than the LO rate due to the 
PDF uncertainties, especially at large m^o. Also, the uncertainty in the total cross section 
at the LHC increases as m^o increases. 

H. Differential cross sections 

Fig. El shows the differential cross section as a function of the transverse momentum pt 
of Z° and A in the associated production of the A°Z° pairs at the LHC. We find that the 
NLO QCD correction could change the shape of transverse momentum distribution. The 
NLO QCD correction enhances the LO differential cross section in low and high px region, 
but reduces in medium p T region. 

Fig. EH shows the differential cross section as a function of the invariant mass M^o^o of 
the A°Z° pairs produced at the LHC. The NLO QCD corrections reduce the LO differential 
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cross sections more in the medium values of M A oz°, and much less in low or high values of 

M A o z o. 

V. CONCLUSIONS 

In conclusion, we have calculated the complete NLO QCD corrections to the inclusive 
total cross sections of the A°Z° pairs produced at the LHC in the MSSM. We have pre- 
formed the calculations using both the DREG and DRED schemes, and found that the 
NLO total cross sections in the above two schemes are the same, which provides a cross 
check to our calculations. Our results show that the LO total cross sections are a few tens 
fb in most of the SUSY parameter space, and can exceed 100 fb for m A o below 160 GeV 
with large tan/3(> 40). The NLO correction can either enhance or reduce the total cross 
sections, but it generally efficiently reduces the dependence of the total cross sections on the 
renormalization/factorization scale. For small m A o and large tan/3, the A"-factor of SUSY 
QCD corrections could become large, and using the QCD plus SUSY improved Yukawa 
coupling in the calculation could reduce the size of the overall if-factor. We have also ex- 
amined the uncertainty in total cross sections due to the PDF uncertainties, and found that 
the uncertainty in NLO cross sections is slightly larger than that in LO ones, especially at 
large m A o. Finally, we also examined a few differential distributions and found that the 
NLO QCD corrections could change the shape of transverse momentum and invariant mass 
distributions. 
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Appendix A 

In this appendix, we give the relevant Feynman rules. 
1. h°(H°) -b-b: A m m b 

A i = n n , M 



2m w cos (3 2m w cos (3 

where a is the mixing angle in the CP even neutral Higgs boson sector |2[. Here, we use the 
abbreviations s a = sin a and c a = cos a. 
2. A - b-b: A 3 m b 'y 5 

-gtaap 
/±3 — . 

2m,„ 



3. h°(H°) -Z°-A°: F m ( Ph0{HO) +p A o] 



gcos([3-a) -gsm((3-a) 

— ^ 7, 5-^2 — 



2 cos 9 W 2 cos 9 W 

Here, we define the outgoing four-momenta of h°(H°) and A to be negative and positive, 
respectively. 

4. z° - b - b : 7„(C V + C Al5 ) 



-w 



1 2 

2 + 3'' 



- 7: + ^sin 2 9 w ), Ca 



-ig 



4 cos 9,, 



2 cos 9,_ 

5. A ) -b a -bp: i[R b G\ %z) {R b ) T U 

j^h b [A b s a + fic a } &(~5 sin ^) sin (« + Z 3 ) + v^mb^s 



V 



V 



/ 



sin 91) cos(« + p) - V2m b h b c a j 



2m 



ir 







Aft tan /3 + /i 




with h b = 

V^mvF cos p 

6. Z° - 6 a - L : 



-|cos^+|sin^ 
±sin2% 



±sin20 s 
-isin^+isin^ 



where and are the four-momenta of b a and bp in direction of the charge flow. 
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Appendix B 



In this appendix, we collect the explicit expressions of the nonzero form factors in Eq. (|14|) 
and Eq. JEJ. Since £M°Mj = 

7 8. ..12 — 0' only the form factors of the first six matrix 
elements are presented here. For simplicity, we introduce the following abbreviations for 
the Passarino-Veltman three-point integrals C^-j and four-point integrals D^j\, which are 
defined similar to Ref. [2^ except that we take internal masses squared as arguments: 



°<Cj) 


= 


(0,0,s,0,0,0), 




rib 
U i{j) 


= 


(m^o, t, 0,0, 0,0), 




nc 


= 


(m|, M ,0,0,0,0), 






= 


(m%,0,u,0, 0,0), 




c e 


= 


K,o,t, 0,0,0), 




u m 


= 


MX, 0,0,0), 




U i{j) 


= 


(u,0, m 2 z , 0,0,0), 




r^h 

U iU) 


= 


{u,0,m 2 AO , 0,0,0), 




fii 
U i(j) 


= 


(m|,t,0,0,0,0), 






= 


( M ,m^ ,0,0,0,0), 




nk 


= 


(m^ o ,t,0,0,0,0), 




nl 

U i(j) 


= 


(0, m|, «, 0,0,0), 






= 


(s,0,0,0,0,0), 






= 


(0,n,m|, 0,0,0), 






= 


(s,m\ ,m%, 0,0,0), 




nv 


= 


{0,m 2 AO ,t, 0,0,0), 






= 


(0,t,m|, 0,0,0), 




nm 


(a, b) = 


C i(i) (0,0,s,mf a ,m|,m 




U i(j) 


(a, b) = 


Ci(j)(m 2 AO ,t,0,ml b ,ml a 




no 


(a, b) = 


C i(j) (0,u,m%,ml b ,ml a 


m|) 


nv 


(a, b) = 


Ci (i )(m^ o ,0,M,m? 6 ,m| 


,> m i 


u m 


(a, b) = 


C i{j) (t, m 2 z , 0, m|, m|, m?J, 




(a, b) = 


Ci (i )(M,0,m^ ,m? 6 ,m| 


,m| 




(a, b) = 


C i( ,)(n,mio,0,m?,m? b 


,m? 


fit 


(a, b) = 


C i(i )(0,m^ o ,t,m?,m? b , 


m\ " 
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D h) = D W ( s > °> *> m A« , 0, m|, 0, 0, 0, 0), 

D iU) = D *U)(°, *> ™l> s > ™ao, 0, 0, 0, 0, 0), 

D i(j) = A(j)(s, 0, u, m|, 0, m^o, 0, 0, 0, 0), 

Df U) = A {J )(0, «, m^o, s, m], 0, 0, 0, 0, 0), 

D i(j)( a , M) = A(i)(s, 0, t, m 2 A o, 0, m%, m\ h ,m\ a ,m\ m? ), 

D i(j)( a > & ' = A(j)(0, 0, m|, m^ , s, f , m\ b , mj, m^, m? ), 

Axj)( a > & > = A(j)( s > °' M ' m l> °' m AQ, m ^ m l a i m % m l)- 

Many of the above functions contain the soft and/or collinear singularities. Since all the 
Passarino-Veltman integrals can be written as a combination of the scalar functions A , B Q , 
Co and D , we present here the explicit expressions for the Co and D functions used in our 
calculations: 



C a o 



C d o 



C c o 



CI 



C b o 



C u o 



C h o 



C 9 o 

c{ 



-Li 



Ce 

s 



7T 

y 



cl 



u 



Co 



m 



AO 



In 



-u 



in 



A" 



In 2 



in 



A" 



In 2 



c. 



1 



/nri nV 

°0 — °0 



Co 



9i 

St 



C f 



t - m 2 A0 



u — rriz L e 

t — m| 
ilnf- 



— 77 



s 

-77 

2^ S 
—77 



7T 

T 



7T 



In 



-t 



mi 



mi 



--ln 2 (-)-- 

2 v -t ; 2 



AO 



1 



m 



AO 



--In 2 

2 



s 



7T" 

Y 



1 i 2 \ n ( mzmA ° ) + 71 " 



(s - m\){m\ - t) 



-t 
+ Li 



st I V s 



l A° 



77 



-Li 



mi 



mi 



Li 



777. 



AO 



A" 



1 



--In 2 

- ty 2 



m 2 z )(m 2 z - 1) 



+ iln 2 
2 

4 



777; 



In 



(m 2 z 



sin 



AO 



t 2 



+ -ln 2 

4 



777^ 

S 
2 



7?), 



AO 



U 



t 



2 



777', 



t 



In 



f sm 2 A0 



In f ^ In ^ 



777; 



A l l2 
+ - n 2 
2 



777 



A° 



/// 



AO 



= t9q = <-> 77, 777 A o <-> 7T7 2 ,), 



where C e = (4vr / u 2 /s)T(l - e)/r(l - 2e). 

For diagrams(a)-(g) in Fig. 2, we get the form factors as following, respectively, 
-Aa s 



fi 



3vt(s- m 2 ho )(s-m 2 Ho 



-{ism^ArF^s - m 2 Ho ) + A 2 F 2 (s - mfj]^ 



+F im - g (s - m 2 Ho )G\{a, tyR^R^C™ + F 2 m g (s - m 2 JG b 2 (a, b^R^C™} 
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Amga s (R b 2 R h b,l ~ ^0,1-^6,2) rp/,i/ ,w 2 \ , 77 Ah/ u\/ 2 M^wi 

oZTI ^ 2 [ u m2 ; [FiG^o, - m g J + F 2 G 2 (o, 6)(s - m ho )]C , 
67i (s — m ho ){s — m Ho ) 

fa 

4 imbM Cy-C A )a S [{i e){2cL + {t + ^ 

~{s-m\- m 2 Ao )C{ 2 + m 2 Ao C f n - B (t, 0, 0)) + m|C* + tC{] 



Am b ga s A 3 T z (a,b)R b al R b bl 



[2c q 00 + uCZ + uC q 22 + (t- m 2 Ao )(C! + C q u + C q + Cf 2 ) 



3ntcos9 w 

+(s - m| - m 2 Ao )(C! + Cf 2 ) + m 2 Ao (Ct + C q u + Cf + Cf 2 ) + *C? 2 ], 
= " 4 ^[ A3 ^ {(1 " e)[2Cf - uC{ 2 + (t- m%)C{ x + m|C? 



-(s-m| 



- m 2 Ao )C{ 2 + m^C^ - fl (f , 0, 0)] + m|Cj + tftf } 

[2 ^ + uCf + uCf 2 + (t - m^ )(C? + C q u + Cf + C| 2 ) 



-2^^-C) {(1 _ e)[2cL + (t _ mX)c[ + ^ 

+(s - m| - m 2 Aa ){C{ - Cl) + m 2 Ao (2C{ - C* 2 ) - B (t, 0, 0)] + m|^} 

2gm b a s A 3 T z (ab)R b al R b bl ^ 
3iTcos6 w t 00 ' 

'^'"■' {n - e)[2Qf + (t - + m 2 z Ct 

+ (s - m 2 z - m 2 Aa ){C{ - C{) + m 2 Ao (2C{ - C|) - Bb(* , 0, 0)] + m|C*} 

3ncos6 w t 00 ' 

e= 4^ 2 a. (£V + Cii)(cf _ c j )> 



+ ^ -4im b A 3 a s + + ^ _ ^ _ u(jh + , _ 2 ^ 

37TM ... 

-m|(C^ - 2Cf) + (s - m| - rr£o)(C? - C*) + m^Cf + B (u, 0, 0)], 



= e- 



'-C V (C*-C*), 



3tt 

3^ K C o - 0, 0)] 
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+e 8im b A 3 a s Cv[uC , + ^ _ ^ _ ^ + _ ^ * 
o7ra 

-m|(C 2 d - 2Cf ) + (s - m 2 z - m 2 A o)(Cf - C%) + m 2 A0 Cf + B (u, 0, 0)], 

~ 2 T^ 3aS (Cv + C ^ m *P" - B (u,0,0)} + ^<^M)i&ite 
+ -2irn b A 3 a s + * + _ ^ _ _ ^ 

OTTU 

-m\C% + m A0 C* + B (u, 0, 0)], 
3^ [™a Co -£ (u,0,0)] 

H LJiL_?_J_i [i?^ 2 i?^ (Cy - Ca)Cq - i^-R^Cy + C A )Cq] 

61XU 

+e ^rn b A 3 a s Cv[uC , + _ _ ? _ ^ 

07TU 

-m|C| + m^oCf + S (u, 0, 0)], 

4i(Cy - C A )m b A 3 a s ^ k _ 2 -4(Cy - CUjn^a^M^A ^ 
— [C Q m Ao - B (t,0,0)] + C 

4i(C v - C A )m b A 3 a s b 2 
3^t [^i m A + B o{t, 0, 0)], 



-8iCvm 6 A3Q! a _pk m 2 A _ B ^ t Q)] + -4(Cy + C A )m~ g a s G b 3 (a, b)R b a2 R b bl ^ 



3vrf ! 3vrf 

37Tt 

A 3 c 

3nt 



A(C A -C v )m- g a s G b (a,b)R b al R b b2 ^ 
_l c ( 

+e ^— [Clm Ao + 5 (t, 0, 0)], 



-n 



2 ' 



2 ' 

-4i(Cy + C A )m b A 3 a s _ [2C , o + ^ + ^ + ^ + ^ 



37TU 



+ ~ 3 4 ^^ T z (a, 6)(2C ° - uC[ 2 ){R b al R b bl + 
Z^V±^^ [2C g o + uCfi + ^ + m | (C c + ^ }] 

+ 3nucose w Tz{amAC ^ 
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i = 2i(Cv ^^ Ri^rnlBoiO, ml ml) + mjBoit, m% m\) - 2m? 



+2mlB (0, m? a , m?J - m\B^ m?, m?J - 2m\ a + tB (t, m?, m\) 
+ 2mlB (0,ml,ml)] + W v - CJm^a. {1 _ ^ M>0)> 

fi = ~ 2i 2p aS Wv ~ C *) R li + (°v + [-2m?S (0, ml ml + m?B (t, mj, m\) 

+2mlB (0, m\ a ,m\ a ) - 2m? - m\B Q {t, m?, m?J - 2m\ a + t£ (t, m?, m?J 

+2mjB (0,mlmi)} + " 2 ^f y ^ (l - e)fBb(f,0,0), 
f! = ~ 2%{CV + 3 ^ )mbA3as Rl 2 2 [-2mlB (0, ml m|) + m^ (t, mj, m?J + 2m?E (0, <, m 



-2m? - ml a B (u, m?, m?J - 2m? a + uS (u, m?, m?J + 2m? aJ B (0, m?, m?J] 

-2i(Cv + C A )m b A 3 a s 

3^2 ~ e ) MjB o(w, 0, 0), 

fi = 2 ^f S [(Cv - ^)< 2 + (Cv + C A )^ 2 2 ][-2m^ (0,m|,m|) + m\B^mlm\ a ) 
+2m~ -BnfO, mr , m\ ) — 2m? — m| -B (u, m? ml ) — 2m? + uB (u, m? m~ ) 

Da v 1 Oa 1 Oa' 9 o a x 9 1 b a ' b a v ' 9 7 b a ' 

+2m? a 5 (0,m|,m?J] + 2im ^C v a s {1 _ c)uSo(li> , 0), 
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f 9 - 

J 5 2 , 

f/ _ _ /2 

; 6 --g- 



For the box diagrams(a)-(d) in Fig. 3, we find, respectively, 
+(l + e)(C£-C*)] 

+e ' W3( 3 C ;' CA)as [C 2 e + Cr + (« - m^)(^ + - (t - m%)(D b 2 + D») 

-rr£o(I>? + + + £&)], 
= 4im ^^-^ Q ' [m! ; nW + D\) - (t - m\ Q ){Dl + D 2 ° + - (u - m\ Q )D\ 
-C a + CI - ci\ 

+e ^MCv^-C A )a S [2D ^ _ {u _ m , o)( ^ + ^ + ^ + ^ 



+mio(^ 3 + ^23 + ^33 + ^2)], 
— 2Cyf B ° X ^ 



Box(a) 

Cy — C/ 

Box(a) 



rBox(a) — 2Cyf 2 
J A ~ " 



Cv — C/ 
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f B 0x{a) = 2im b MCv - C A )a, + _ + + _ 

37T 

+e 2^ 3 (Cv--C) a , [20gn _ (i _ _ m%DI _ ^ 

nn f Box(a) 
rBox(a) _ ~ ih 

^ = JZZU. ™ a)6 ^ - w + ^ + 

= -4,:m t/t3 (C,, + C, )a , | _ sg; + (u _ m|)(fl , _ Dl) _ m%D} _ c;] 
+e -W 3 (CV + C,)a, [c? _ c „ _ Cf _ (( _ m|)fl;i + (u _ m%)Dl 
+m%(D\ + £>«,)], 

f B„ (c) = -4,m„A 3 (Cv + C A ) a , [cr _ _ _ {t _ ml)Di _ {u _ m%)w + ^ 
+m%(D'„ + Dt) - C$] 

+e -W 3 (CV + C,)a, [cf + c , ^ (f m|)(D; _ + ^ + ^ 
+mi(£>; 3 + £>| 3 + Bf) + CJ], 

or y f Box(b) 
Box(c) _ —^VJl 



f: 



3 



4 



Box(c) 
6 



-2Cvf B ° x ^ 



Cv + Ca 

/5 ^ (c) = 2»m 6 A 3 (Cv + g i 4)a a ^ + ^ _ 

37T 



+e 2, mt .4 3 (C, + C,) a , |2g;o _ (u _ m2<o)B , _ ^ _ 0J]i 



-2C V /, 



Box(b) 



Cy + Ca 



/f"*" = f^r z (i, tt )Gl(i,&)(<ii - i?;X)(r>{ + b5 
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FIG. 1: Leading order Feynman diagrams for bb — ► A°Z°. 
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FIG. 2: One-loop virtual diagrams, including self-energy and vertex corrections for bb — > A°Z°. 
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FIG. 3: Box diagrams for bb -► A°Z°. 




FIG. 4: Feynman diagrams for the real gluon emission contributions. 
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FIG. 5: Feynman diagrams for the emission of a massless bottom quark contributions. 
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FIG. 6: LO total cross sections of pp — ► A°Z° via bb annihilation, compared with the ones from 
gluon fusion and Drell-Yan processes at the LHC, as a function of tan j3 with m&(m&) = 4.25 GeV, 
assuming: (1) m = 200 GeV, mi = 160 GeV, A = 100 GeV and /z < 0; (2) m = 150 GeV, 
mi = 180 GeV, A = 300 GeV and p > 0. 
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FIG. 7: LO total cross sections of pp — > A°Z° via 66 annihilation, compared with the ones from 
gluon fusion and Drell-Yan process at the LHC, as a function of m A o with m&(mb) = 4.25 GeV, 
assuming: mi = 160 GeV, A = 100 GeV, tan/3 = 40 and \x < 0. 
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FIG. 8: Dependence of the NLO total cross sections for the A°Z° production at the LHC on 
the theoretical cutoff cale 6 S with 6 C = 5 s /50, assuming: (1) mo = 200 GeV, mi = 160 GeV, 
A = 100 GeV, tan/? = 40 and fi < 0; (2) m = 150 GeV, mi = 180 GeV, A = 300 GeV, 

2 

tan/3 = 40 and \x > 0. Here, we take mb(rrib) = 4.25 GeV. In (a), the solid and dotted curves are 
the results for model (1) and (2), respectively. 
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FIG. 9: Dependence of the total cross section of the A°Z° production at the LHC on m A o, assuming 
mi = 160 GeV, A = 100 GeV, and n < for tan (3 = 40 in Fig. 9(1) and tan = 10 in Fig. 9(2). 
Three differet calculations were done by using: (I) MS bottom quark mass at the scale m^, (II) 
QCD improved bottom quark running mass at the scale m A o, and (III) QCD plus SUSY improved 
bottom quark running mass at the scale m A o, respectively, to evaluate the bottom quark Yukawa 
coupling. 
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FIG. 10: i-C-factor, defined as (Jnlo/vlo, for the A°Z° production at the LHC as a function of m A o, 
using the QCD improved running to evaluate the bottom quark Yukawa coupling, assuming 
mi = 160 GeV, A = 100 GeV, and n < for tan/3 = 40 in Fig. 10(1) and tan/3 = 10 in Fig. 10(2). 

2 

The full K- factor is shown as curve (a), which includes the pure QCD corrections, shown as curve 
(b), and SUSY QCD corrections, shown as curve (c). The contribution from the SUSY QCD box 
diagrams is also separately shown as curve (d) for comparison. 
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FIG. 11: K- factor, defined as <jnlo/&lo, for the A°Z° production at the LHC as a function of m^o, 
using the QCD plus SUSY improved bottom quark Yukawa coupling, assuming mi = 160 GeV, 

2 

Aq = 100 GeV, /i < and tan/3 = 40. The full K- factor is shown as curve (a), which includes 
the pure QCD corrections, shown as curve (b), and SUSY QCD corrections, shown as curve (c). 
The contribution from the SUSY QCD box diagrams is also separately shown as curve (d) for 
comparison. 
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FIG. 12: Comparison of the SUSY QCD corrections, denoted as SasusYQCD/o'LO, f° r the A°Z° 
production at the LHC. The results of using the complete numerical calculation (dashed curves) 
and the approximate analytical forms (solid curves) in the heacvy mass limit are separately shown 
as a function of Msusy with tan/3 = 4 and 40, respectively, assuming m A o = 150 GeV and Mq = 
Mp = fi = Ab = Mg = Msusy- Here, the LO cross section is calculated by using the MS bottom 
quark mass. 
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FIG. 13: Dependence of the total cross sections for the A°Z° production at the LHC on tan/3, 
assuming: (l)m = 200 GeV and 400 GeV, respectively, mi = 160 GeV, A = 100 GeV, and fj, < 0; 
(2)m = 150 GeV, mi = 180 GeV and 400 GeV, respectively, A = 300 GeV, and ft > 0. 
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FIG. 14: Dependence of the if-factor, denned as cfnlo/&lo, on tan/3 for the A°Z° production at 
the LHC, assuming: (l)m = 200 GeV and 400 GeV, respectively, mi = 160 GeV, A = 100 GeV, 
and n < 0; (2)m = 150 GeV, mi = 180 GeV and 400 GeV, respectively, A = 300 GeV, and fi > 0. 
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FIG. 15: Dependence of the total cross sections on the renormalization/factorization scale (fi r = 
fj, f ) for the A°Z° production at the LHC, assuming: (1) = 200 GeV, mi = 160 GeV, A = 

2 

100 GeV, tan/3 = 40 and /z < 0; (2) m = 150 GeV, mi = 180 GeV, A = 300 GeV, tan/3 = 40 

2 

and (i > 0. Here, the QCD plus SUSY improved bottom quark Yukawa coupling is used. The case 
of the curves (3) is similar to (1), but in (3) the pure QCD running bottom quark mass is used 
instead. The case of the curves (4) is similar to (1), but in (4) the contribution from the SUSY-EW 
correction in the running bottom quark Yukawa coupling is not included, namely, only the pure 
QCD and SUYSY-QCD corrections are included. 
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FIG. 16: Dependence of the total cross sections on the factorization scale (///), labelled as case (1), 
or renormalization scale (/U r ), labelled as case (2), for the A°Z° production at the LHC, assuming: 
m = 200 GeV, mi = 160 GeV, A = 100 GeV, tan/3 = 40 and /z < 0. Here, the QCD plus SUSY 

2 

improved bottom quark Yukawa coupling is used and m av = (m^o + m z o)/2. 
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FIG. 17: The PDF dependence of the total cross sections for pp — > A°Z° production at the LHC, 
as a function of Aq, assuming tuq = 250 GeV, mi = 160 GeV, tan/? = 40 and \i < 0. Here, the 

2 

QCD running bottom quark mass is used to evaluate the bottom quark Yukawa coupling. 
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FIG. 18: The PDF dependence of the total cross sections for pp — ► A°Z° at the LHC as a function 
of m^o, assuming Ao = 100 GeV, mi = 160 GeV, tan/3 = 40 and fi < 0. Here, the QCD running 

2 

bottom quark mass is used to evaluate the bottom quark Yukawa coupling. 
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FIG. 19: Differential cross sections in the transverse momentum (px) of Z° and A bosons, for 
the A°Z° production at the LHC, assuming: m = 200 GeV, mi = 160 GeV, A = 100 GeV, 

2 

tan (3 = 40 and \i < 0. 
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FIG. 20: Differential cross sections in the invariant mass (M^o^o), for the A°Z° production at the 
LHC, assuming: (1) m = 200 GeV, mi = 160 GeV, A = 100 GeV, tan/? = 40 and p < 0; (2) 
m = 150 GeV, mi = 180 GeV, A = 300 GeV, tan (3 = 40 and n > 0. 
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